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Abstract.  The  purpose  of  the  present  investigation  is  to 
study  possible  electromagnetic  effects  in  the  Earth’s 
ionosphere  and  magnetosphere,  which  are  caused  by  a 
cosmic  body  impact  onto  the  atmosphere  and  by  a 
cosmic  body  passage  through  the  magnetosphere.  An 
analysis  of  a  cosmic  body  impact  onto  the  atmosphere 
is  carried  out  for  a  meteoroid  with  the  kinetic  energy 
of  the  order  of  1  kt  TNT.  Such  meteoroids  are  usually 
detected  in  the  Earth’s  atmosphere  several  times  a  year. 
Variation  of  the  electron  density  in  the  ionospheric 
plasma  and  excitation  of  plasma  turbulence  as  a  result 
of  meteoroid  impact  events  are  shown  to  be  associated 
with  an  appearance  of  the  acoustic  perturbation  for¬ 
med  in  the  atmosphere  due  to  the  impact  of  the  meteor¬ 
oid.  The  evolution  of  the  acoustic  perturbation  is  stud¬ 
ied.  The  electron  density  variation  and  parameters  of 
ion-acoustic  waves  (and,  in  particular,  ion-acoustic 
sohtons),  which  are  associated  with  the  acoustic  per¬ 
turbation,  are  estimated.  Variations  of  geomagnetic 
field  strength  caused  by  the  meteoroid  impact  events 
and  the  rates  of  electron  acceleration  as  a  result  of 
the  interaction  of  the  cosmic  body  with  the  Earth’s 
atmosphere  are  evaluated.  Most  significant  elec¬ 
tromagnetic  effects  in  the  Earth’s  magnetosphere  can 
be  caused  by  a  cosmic  body  embedded  in  a  dust  cloud. 
It  is  demonstrated  that  for  an  adequate  description  of 
the  dust  cloud  with  the  magnetosphere  the  process  of 
dust  particle  charging  during  the  meteoroid  passage 
should  be  taken  into  account.  The  presence  of  charged 
dust  particles  can  result  in  the  formation  of  a  current 
system  and  (for  sufficiently  large  density  of  dust 
particles,  e.g.  close  to  that  related  to  the  comet  Shoe¬ 
maker-Levy  9)  in  violation  of  the  magnetic  structure 
of  the  magnetosphere  in  the  vicinity  of  the  dust  cloud. 
The  possibility  of  auroral  processes  before  the  effective 
interaction  of  the  cosmic  body  with  the  atmosphere, 
which  are  triggered  by  fast  particles  created  due  to 
Fermi  acceleration  process,  is  demonstrated.  ©  1997 
Elsevier  Science  Ltd 


Introduction 

Some  large  meteoroids  detected  by  satellites  in  the  Earth's 
atmosphere  (Tagliaferri  et  al.,  1994)  have  a  kinetic  energy 
of  the  order  of  1  kt  TNT1  (Adushkin  and  Nemchinov, 
1994).  Such  meteoroids  are  usually  detected  several  times 
a  year.  In  this  article  we  shall  concentrate  most  attention 
on  such  meteoroids.  We  note,  however,  that,  although 
larger  meteoroids  (with  kinetic  energy  of  the  order  of  1  Mt 
TNT)  occur  more  rarely  in  the  atmosphere,  their  detection 
is  possible  rather  often  (one  event  per  10-20yr).  Some 
satellite  measurements  and  those  carried  out  on  the  Earth 
have  shown  the  possibility  of  the  appearance  of  oscil¬ 
lations  in  the  magnetospheric  and  ionospheric  plasmas 
over  earthquake  centres,  explosive  volcano  eruptions, 
explosions  in  the  atmos'phere,  etc  (Gokhberg  et  al.,  1 982 ; 
Gohberg  et  al.,  1982;  Sorokin  and  Fedorovitch,  1982; 
Galperin  et  al.,  1985).  Plasma  turbulence  generated  as  a 
result  of  perturbations  in  the  atmosphere  (e.g.  due  to 
explosions)  can  be  related  to  a  change  in  the  electron 
density  variations  in  the  ionosphere.  This  is  confirmed,  in 
particular,  by  analysis  (Nagorsky,  1985)  of  consequences 
of  the  explosion  carried  out  within  the  frames  of  the  pro¬ 
ject  MASSA  (with  the  energy  release  of  260 1  TNT).  Thus 
one  can  expect  (analogously  to  the  case  of  usual 
explosions)  that  the  energy  release  in  the  process  of  an 
impact  of  a  meteoroid  with  high  kinetic  energy  is 
accompanied  by  a  change  in  the  electron  density  in  the 
ionosphere.  This  effect  is,  in  particular,  caused  by  the 
influence  on  the  ionospheric  plasma  of  an  acoustic  per- 


1  In  the  present  paper  we  use  the  quantity  kt  TNT  (kilotons 
TNT)  to  denote  the  kinetic  energy  of  the  meteoroid.  1  kt  TNT 
corresponds  to  4.2  x  10l')  erg.  This  quantity  is  usually  used  in  the 
context  of  nuclear  explosions  and  is  often  used  in  the  description 
of  meteoroids  (Adushkin  and  Nemchinov,  1994).  A  kinetic- 
energy  of  1  kt  TNT  for  a  meteoroid  moving  with  speed 
2xl0<’cms  _l  corresponds  to  a  mass  of  2.1  v  10  g  of  the 
meteoroid. 
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turbation  formed  in  the  process  of  the  energy  release.  The 
electron  density  variations  in  the  ionosphere  can  sig¬ 
nificantly  influence  radio  wave  propagation  determining 
violation  of  phases,  azimuth  angles  of  the  reflected  waves, 
the  characteristics  of  multiple-beam  radio  wave  propa¬ 
gation,  etc.  Furthermore,  the  electron  density  variations 
and  acoustic  perturbation  caused  by  the  meteoroid  impact 
can  result  in  excitation  of  plasma  turbulence  in  the  iono¬ 
sphere.  Generation  of  plasma  turbulence  can  lead  to 
sufficiently  effective  electron  acceleration  due  to  their 
interaction  with  waves  in  a  plasma.  This,  in  turn,  can 
result  in  excitation  of  atoms  and  molecules  at  altitudes 
exceeding  90  km  in  the  region  of  the  ionosphere  mag¬ 
netically  related  to  the  region  of  the  energy  release  and, 
as  a  consequence,  in  glow  of  this  region.  If  a  cosmic 
body  is  embedded  in  a  dust  cloud  (such  as  the  comet 
Shoemaker-Levy  9  (Ip  and  Prange,  1994;  Dulk  and 
Leblanc,  1995))  and  moves  through  the  magnetosphere  of 
the  Earth  then  the  presence  of  the  dust  can  also  lead  to 
electromagnetic  effects  in  the  ionosphere  and  mag¬ 
netosphere. 

Here,  we  estimate  the  effect  of  a  change  in  the  electron 
density  in  the  ionosphere  in  the  process  of  propagation  of 
the  acoustic  perturbation  and  discuss  possible  difference 
between  its  manifestations  for  usual  explosions  and  for 
meteoroid  events.  We  discuss  the  possibility  of  excitation 
of  ion-acoustic  oscillations  in  the  ionospheric  plasma  by 
the  acoustic  perturbation  and  possibility  of  an  appearance 
of  ion-acoustic  solitons  in  the  ionospheric  plasma.  We 
present  the  main  parameters  characterizing  the  ion-acous¬ 
tic  turbulence  caused  by  the  meteoroid  events.  We  men¬ 
tion  briefly  the  possible  influence  of  the  meteoroid  events 
on  the  geomagnetic  activity  and  give  some  estimates  of 
the  rates  of  electron  acceleration  as  a  result  of  the  inter¬ 
action  of  the  cosmic  body  with  the  Earth’s  atmosphere. 
Finally,  we  consider  effects  caused  by  the  presence  of  dust 
particles  if  the  cosmic  body  moving  through  the  Earth’s 
magnetosphere  is  embedded  in  a  dust  cloud. 


Acoustic  wave 

Variation  of  the  electron  density  in  the  ionospheric  plasma 
and  excitation  of  plasma  turbulence  as  a  result  of  the 
meteoroid  events  are  connected  with  an  appearance  of  an 
acoustic  perturbation  in  the  atmosphere.  We  assume  that 
the  evolution  of  the  acoustic  perturbation  formed  in  the 
atmosphere  due  to  the  impact  of  a  meteoroid  with  the 
kinetic  energy  of  the  order  of  1  kt  TNT  is  the  following : 
the  acoustic  shock  perturbation  is  formed  at  the  distance 
r0  from  the  region  of  maximum  energy  release ;  the  form 
of  the  perturbation  profile  can  be  approximated  initially 
as  a  Glasstone  pulse  (Glasstone  and  Dolan,  1977)  ;  later 
the  rarefied  wave  of  this  pulse  is  transformed  to  the  second 
shock  wave  and  the  Glasstone  pulse  is  converted  to  an  N- 
like  wave.  Such  evolution  is  analogous  to  that  described 
by  Nagorsky  (1985)  for  the  case  of  usual  TNT  explosions. 

In  analysing  the  evolution  of  the  N-like  wave,  we  sup¬ 
pose  that  the  averaged  wind  is  absent  and  neglect  the 
temperature  gradient  as  compared  with  pressure  and  den¬ 
sity  gradients,  but,  we  take  into  account  the  dependence 
of  the  parameters  of  the  problem  on  the  temperature. 


Furthermore,  for  simplicity  we  use  the  approximation 
of  short  acoustic  waves  (Romanova,  1970)  and  assume 
Knudsen  and  Mach  numbers  to  be  small.  Under  these 
assumptions  the  evolution  of  the  N-like  wave  is  described 
by  the  Navier-Stokes  set  of  equations  where  the  gravi¬ 
tational  effects  should  be  taken  into  account,  the  equa¬ 
tions  of  state,  and  statics  equations.  The  expansion  of  the 
unknown  functions  in  terms  of  Mach  number  in  the  above 
equations  and  consideration  in  the  frame  of  reference 
moving  together  with  the  perturbation  allow  us  to  derive 
an  evolution  equation  for  the  N-like  wave  which  is  anal¬ 
ogous  to  Burgers  equation  (Nagorsky,  1985) : 


y  =  cpjcv,  n  =  0,  1 ,  2  characterizes  the  extent  of  per¬ 
turbation  divergence,  r  is  the  radial  distance  from  the 
explosion  centre,  V0  is  the  ratio  of  the  amplitude  of  the 
oscillation  speed  of  gas  particles  in  the  shock  wave  to  the 
acoustic  speed  c  at  r  —  r0,  V  =  V(r)  is  the  compressional 
pulse  profile,  n 2(r)  is  the  length  of  the  compressional  pulse, 
X,  =  A(r0),  cp  is  the  zenith  angle,  v0  is  the  damping 
coefficient,  r  is  the  phase  coordinate,  H  =  k  T/Mag  is  the 
height  of  the  homogeneous  atmosphere,  k  is  the  Boltz¬ 
mann  constant,  T  is  the  temperature,  M„  is  the  atomic 
mass  of  the  neutral,  g  is  the  gravity  acceleration. 

The  solution  to  equation  (1)  shows  that  the  evolution 
of  the  acoustic  perturbation  depends  significantly  on  its 
geometry.  The  usual  TNT  explosions  result  in  a  three- 
dimensional  “cupola-shaped”  form  of  the  acoustic  per¬ 
turbations  (n  =  2  in  equation  (1)).  The  geometry  of  the 
acoustic  perturbation  caused  by  the  energy  release  in  the 
process  of  the  meteoroid  impact  is  more  complicated.  This 
is  because  the  meteoroid  impact  onto  the  atmosphere 
cannot  be  considered  as  a  spheric  explosion  with  a  definite 
explosion  centre.  At  least,  at  the  first  stage  this  process 
has  cylindrical  geometry  (n  =  1  in  equation  (1)).  Unfor¬ 
tunately,  the  problem  of  formation  and  development  of 
acoustic  perturbation  caused  due  to  the  meteoroid  impact 
onto  the  atmosphere  is  not  still  completely  solved.  In 
particular,  it  is  still  not  clear  at  which  altitudes  the  form 
of  the  acoustic  perturbation  loses  the  cylindrical  features 
and  the  perturbation  can  be  considered  as  three-dimen¬ 
sional  “cupola-shaped”.  This  is  connected  with  the  fact 
that  the  characteristic  size  of  the  region  of  the  meteoroid 
impact  is  the  size  of  the  homogeneous  atmosphere,  thus 
the  cylindrical  features  in  the  form  of  the  acoustic  per¬ 
turbations  can  be  conserved  over  rather  long  period  of 
time.  We  emphasize  that  an  analysis  (for  cylindrical 
geometry  of  the  energy  release)  of  the  generation  and 
propagation  characteristics  of  infrasonic  pressure  waves 
excited  during  meteoroid  entry  into  the  Earth’s  atmo¬ 
sphere  has  been  carried  out  by  ReVelle  (1976).  Unfor¬ 
tunately,  the  final  solution  to  the  problem  of  whether 
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the  cylindrical  features  of  the  acoustic  perturbation  are 
conserved  at  the  altitudes  which  are  of  interest  from  the 
Viewpoint  of  study  of  electromagnetic  effects  has  not  yet 
fseen  proposed.  This  problem  requires  further  inves¬ 
tigations. 

In  the  following  we  use  for  estimates  of  the  charac¬ 
teristics  of  the  acoustic  perturbation  and  of  a  change  in 
electron  density  caused  by  this  perturbation  the  results 
obtained  under  the  assumption  of  a  three-dimensional 
“cupola-shaped”  form  of  the  perturbation.  This  is  related 
to  the  fact  that  we  are  interested  in  altitudes  A  — 120— 
200  km  which  significantly  exceed  the  height  H  of  the 
homogeneous  atmosphere.  One  can  expect  that  the  influ¬ 
ence  of  the  cylindrical  features  in  the  geometry  of  the 
perturbation  is  not  significant  at  these  altitudes.  Fur¬ 
thermore,  estimates  fulfilled  for  the  cylindrical  geometry 
of  the  acoustic  perturbation  show  us  that  the  cylindrical 
features  in  the  geometry  of  the  perturbation  result  in  a 
weaker  acoustic  perturbation  than  the  “cupola-shaped” 
one  due  to  the  usual  explosion  with  the  same  energy 
release.  Thus  the  use  of  the  model  of  “cupola-shaped” 
geometry  (generally  speaking,  with  a  reduced  central 
energy)  is  legitimate  for  the  description  of  the  acoustic 
perturbation  in  our  case,  when  we  consider  the  altitudes 
h  ~  120-200  which  are  of  interest  from  the  viewpoint  of 
the  study  of  electromagnetic  effects. 

For  energy  releases  of  the  order  of  1  kt  TNT  we  assume 
that  r0~105cm,  the  ratio  of  the  amplitude  of  the  oscil¬ 
lation  speed  of  gas  particles  in  the  shock  wave  to  the 
acoustic  speed  c  at  r  =  r0  is  F0~5  x  10-2,  and  the  length 
of  the  compressional  pulse  (at  r  =  r0)  is  nk0  ~  1 04  cm.  These 
parameters  are  close  to  those  used  by  Nagorsky  (1985)  to 
describe  the  usual  explosion  of  energy  release  of  260 1 
TNT  carried  out  on  the  Earth.  Furthermore,  these  data 
allow  us  to  find  the  values  of  the  relative  amplitude  (A) 
and  the  length  (nX)  of  the  compressional  pulse  at  the 
altitudes  of  120-1 50  km  (where  the  maximum  values  of 
the  wave  amplitude  are  attained)  which  are  of  the  order 
of  the  magnitudes  of  A  and  X  obtained  (for  the  energy 
release  of  1  kt  TNT  occurring  at  h0  ss  40  km)  on  the  basis 
of  formulae  (Gubkin,  1958,  1961)  which  describe  the 
dependence  of  the  relative  amplitude  on  the  distance  R 
from  the  energy  release  centre  and  on  the  energy  released. 
Here  we  present  Gubkin’s  formula  for  the  case  of  an 
acoustic  perturbation  of  three-dimensional  (“cupola- 
shaped”)  geometry  and  for  vertical  propagation  of  the 
perturbation : 

Ap  0.135£"/3  \hQ  R) 

A~  P  ~  J?72-68  +  °-17*-ln£,/3eXP  13//  +  2Hi 

where  E  is  the  energy  released  in  the  meteoroid  impact 
event  expressed  in  units  of  kt  TNT ;  R,  ho,  and  H  are 
calculated  in  units  of  km ;  p  is  the  pressure  at  the  altitude 
h0  +  R ;  and  Ap  is  the  pressure  variation  in  the  wave.  The 
presence  of  the  maximum  values  of  the  wave  amplitude 
at  altitudes  of  about  120-1 50  km  is  connected  with  the 
fact  that  the  main  energy  absorption  of  the  acoustic  per¬ 
turbation  is  related  to  the  non-linear  absorption  for  the 
altitudes  below  120-1 50  km.  At  the  altitudes  of  120— 
180  km  the  decrease  in  perturbation  amplitude  is  also 
related  to  the  increase  of  the  acoustic  velocity  with  the 
altitude.  Above  170-1 80  km  the  main  contribution  to  the 


amplitude  decrease  is  made  by  the  viscosity.  The  length 
of  the  compressional  pulse  increases  with  altitude  so  that 
at  an  altitude  of  about  200  km  the  length  of  the  pulse  is 
two  orders  of  magnitude  larger  than  at  the  surface  r  =  r0. 
For r0~  105cm,  F0~5  x  10-2,  nX0~  104cm, and hn x 40 km 
we  find  that  at  the  altitude  h ~  150 km  A~2x  10_l. 
nX~  5  x  105cm. 

We  emphasize  that  the  above  results  can  be  influenced 
by  the  presence  of  the  wake  behind  the  meteoroid.  The 
gas  in  this  wake  has  the  characteristics  different  from  those 
of  the  ambient  gas.  Accordingly,  the  properties  of  the 
acoustic  perturbation  propagating  inside  the  wake  are 
different  from  those  in  the  ambient  gas. 


Electron  density  variation 

Estimate  the  electron  density  perturbation  in  the  iono¬ 
sphere  caused  by  the  meteoroid  events.  We  present  the 
electron  density  in  the  form  ne(r,t)  =  n0(h)[\+n'e(r,t)], 
where  n0  is  an  unperturbed  value  of  the  electron  density, 
n'c(r,t)  is  the  electron  density  deviation  from  the  unper¬ 
turbed  state  caused  by  the  acoustic  perturbation.  We  shall 
consider  separately  three  regions : 

•  the  region  E  of  the  ionosphere  where  <wBi «  vi„ ; 

•  the  lower  part  of  the  region  EF  where 
n0~  const; 

•  the  upper  part  of  the  region  EF  and  the  region  F  of  the 
ionosphere  where  a>Bi»  vin. 

Here  a>Bj  is  the  ion  gyrofrequency,  vin  is  the  frequency  of 
binary  collisions  between  ions  and  neutral  particles. 

Under  the  assumption  that  the  electron  density  dis¬ 
tribution  in  the  region  E  is  given  by  the  expression 
n0««mexp[— (h— /jm)2/AA2]  (where  nm  is  the  maximum 
value  of  the  electron  density  in  the  region  E,  Ah  is  the 
characteristic  width  of  this  region)  the  electron  density 
variation  in  the  region  E  is  (Gokhberg  et  al.,  1982) 

n'  =  V—2  A  —j^r  W  cos  <P  (2) 

where 

'(t-t2/2)/tt  ifO<T<7t 
A  In  cosh(u/ A)  +  x  if  t  ~  0 . 

In  the  lower  part  of  the  region  EF  we  have  (Gershman, 
1974) 

K  =  y—^rf  1  +  ^cos2  (4) 

«B,  +  vi\  V2n  ) 

where  9  is  the  angle  between  the  wave  vector  k  of  the 
acoustic  perturbation  and  the  direction  of  the  geo¬ 
magnetic  field. 

In  the  region  F  of  the  ionosphere  the  electron  density 
perturbation  in  the  frame  of  reference  created  by  the  axes 
Ox  1,  Oyu  Ozi  with  the  axis  Ozx  in  the  direction  parallel 
to  the  geomagnetic  field  is  described  by  the  expression 
(Nagorsky,  1985) 
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\-J^rai2"zi,'~nc l,'dz;  <5) 

where  G(zuz[,t— t')  is  the  Green  function  given  by 

G(zx,z[,t-t')  =  exp  (6) 

the  subscript  ||  denotes  the  vector  component  parallel  to 
the  geomagnetic  field,  D  =  k(Tc+  T)IM\v-m  is  the 
coefficient  of  the  ambipolar  diffusion,  Mx  is  the  ion  mass, 
re(i)  is  the  electron  (ion)  temperature.  Equation  (5)  is 
obtained  under  the  assumption  that  the  ionosphere  is  a 
layer  with  the  width  L  (so  that  L »  2n/k^  and  that  the 
boundaries  of  the  layer  do  not  influence  the  deviation  of 
the  electron  density  from  the  steady-state  values.  This 
enables  to  set  the  boundary  condition  n'e  =  0  at  infinity. 
The  Fourier  transformation  for  the  values  d,  and  «'  gives 
the  following  result  for  the  Fourier  component  of  the 
electron  density  variation 


where  a  =  (n02 — tto^n/rtok^L,  n0H2)  is  the  unperturbed  value 
of  the  electron  density  at  the  lower  (upper)  boundary  of 
the  ionospheric  layer,  q  =  k2K(re+  T^lMfijvm. 

Thus  the  electron  density  variation  is  proportional  to 
the  amplitude  A  of  the  acoustic  perturbation  in  all  three 
regions  considered  of  the  ionosphere.  For  the  meteoroid 
impact  resulting  in  the  acoustic  perturbation  which  is 
characterized  by  n  =  2,  r(j~  105 cm,  K0~5  x  10~2,  and 
nkQ~  104cm  the  maximum  values  of  the  electron  density 
variations  n'e  are  attained  at  the  altitudes  h  » 120-1 50  km 
and  are  of  the  order  of  0.1.  If  the  cylindrical  features  in 
the  geometry  of  the  acoustic  perturbation  are  significant 
up  to  the  altitudes  h  ~  1 50-200  km  (so  that  we  can  use  the 
estimate  n  =  1,  r0~  102-103  cm,  jtlo<102cm,  and 

V0  ~  5  x  10-2)  then  the  maximum  value  of  the  variation  of 
the  electron  density  is  smaller :  «' «  0. 1 . 

We  emphasize  that  the  above  estimates  of  the  electron 
density  variation  imply  that  the  regions  E  and  F  of  the 
ionosphere  are  not  stratified.  However,  an  appearance  of 
different  layers  in  these  regions  can  lead  to  significant 
changes  in  both  the  form  and  the  amplitude  of  the  acoustic 
perturbation.  Furthermore,  the  model  used  implies  that 
the  acoustic  perturbation  propagates  in  an  atmosphere 
gas  consisting  of  neutral  particles.  However,  in  the  iono¬ 
sphere  the  usual  acoustic  perturbation  can  result  in  the 
generation  of  ion-acoustic  waves  and  other  plasma  modes 
and  through  this  it  can  be  modified.  Moreover  the  plasma 
modes  can,  in  turn,  influence  the  acoustic  perturbation. 
All  these  effects  can  modify  the  above  estimates. 


Oscillations  in  plasma 

We  now  discuss  the  possible  appearance  of  oscillations  in 
the  ionospheric  plasma  due  to  the  effect  of  the  acoustic 


perturbation  created  in  the  process  of  impact  of  the 
meteoroid.  Here  we  restrict  ourselves  to  the  case  when  the 
effect  of  the  energy  release  in  this  process  can  be  described 
in  terms  of  an  explosion  with  n  =  2,  r0~105cm, 
F„~5x  10-2,  and  n/.~  104cm,  i.e.  when  we  can  neglect 
the  cylindrical  features  in  the  form  of  the  acoustic  per¬ 
turbation. 

As  has  been  shown  above  the  acoustic  perturbation 
created  in  the  process  of  impact  of  the  meteoroid  can 
result  in  an  appearance  of  density  variations  n'  ~  0. 1 .  This 
effect  is  connected  with  the  motion  of  charged  particles 
arising  due  to  collisions  of  neutral  particles  moved  by 
the  acoustic  perturbation  with  the  charged  particles.  The 
acoustic  perturbation  can,  in  particular,  be  converted  to 
the  ion-acoustic  mode.  This  occurs  at  altitudes  of  1 50— 
200  km  (the  so-called  altitudes  of  the  acoustic  zone)  where 
the  frequency  of  inelastic  collisions  exceeds  that  of  elastic 
collisions  (Alebastrov  et  at.,  1985).  The  time  at  which  the 
acoustic  perturbation  reaches  these  altitudes  is  of  the 
order  of  5  x  102s.  Since  in  the  process  of  the  mode  con¬ 
version  the  frequency  co  is  assumed  to  be  constant,  we  can 
estimate  the  frequency  of  the  ion-acoustic  mode  and  the 
component  parallel  to  the  geomagnetic  field  of  its  wave 
vector.  Indeed,  the  characteristic  value  of  the  length 
of  the  acoustic  pulse  at  the  altitude  /i  ~  1 50  km 
is  5  x  105cm,  the  acoustic  velocity  is 

c  =  JyKTjMa~5x  1 04 cm s ~ 1 ,  the  ion-sound  velocity  is 
vis  =  JkTM ~  5  x  104cms_1  (we  suppose  that  y  =  5/3, 
r ~  600  K ;  the  electron  temperature  Tt  ~  800  K ;  the  con¬ 
stituents  of  the  plasma  at  this  altitude  are  electrons,  the 
ions  NO+  and  02+,  and  neutral  particles,  the  main  neutral 
components  are  N2, 0,  and  02  (Mishin  et  al.,  1989)).  Thus 
the  characteristic  value  of  the  wave  vector  length  of  the 
acoustic  perturbation  is  &a~ 2n/k *4  x  10-5 cm- 1  and  the 
characteristic  frequency  is  co  =  £ac~  1  s- '.  This  frequency 
is  significantly  less  than  the  ion  gyrofrequency  at  the  alti¬ 
tudes  considered:  coBi  — 102 s' '.  Furthermore,  the  fre¬ 
quency  co  is  much  less  than  the  frequency  of  collisions 
between  electrons  and  other  plasma  particles.  The  latter 
is  ve ~  ven ~  vin ~  103  s~ 1  at  h~  150  km  (Mishin  et  al.,  1989). 
Thus  the  acoustic  perturbation  is  converted  to  the  mag¬ 
netized  ion-acoustic  mode  in  a  collisional  plasma,  i.e.  its 
dispersion  law  is  co  =  ktcs,  where  cs  is  the  ion  sound  speed 
for  such  a  plasma  which  differs,  in  general,  from  vis  (see, 
e.g.  Vladimirov  et  al.,  1995).  We  note  that 

M»3^ve  (8) 

where  Mt  is  the  electron  mass.  The  characteristic  fre¬ 
quency  3veMJM-,  corresponds  to  the  inverse  time  of  the 
electron  and  ion  temperature  equalization.  Consequently, 
we  can  assume  that  in  the  wave  the  relationship  between 
the  electron  (Tc)  and  ion  (7))  temperatures  is  conserved, 
i.e.  since  at  /i~150km  Tz  ~  7j,  we  can  consider  that  the 
latter  relation  remains  valid  in  the  wave.  In  this  case  we 
have  (Vladimirov  et  al.,  1995,  p.  334) 


r  Me/  5rAT/2 


where  t>Te  =  kTJMc  is  the  electron  thermal  velocity. 
Thus  the  component  parallel  to  the  geomagnetic  field  of 
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the  wave  vector  kis  of  the  ion-acoustic  wave  is  of  the  order 
of  the  characteristic  length  of  the  wave  vector  ka  of  the 
acoustic  perturbation :  |A:iS|||~|ka|ss4  x  10~5cm_1. 

The  collisional  damping  rate  for  these  waves  is  (Vla¬ 
dimirov  et  al .,  1995)  ys  —  y(e)+ y(,),  where 


1.92  +  0.64(7;/'ri)  k\v\ s 
(5/3 )  +  (TJTd  v,  ' 


(10) 


From  these  formulae  we  can  estimate  the  characteristic 
time  of  existence  of  the  ion-acoustic  wave:  Tis~|ysr‘. 
Assuming  that  M,=  M0+,  vc~  103s_l,  and  estimating 
ys  =  y<e),  we  find  tis  ~  200  s.  This  time  is  significantly  longer 
than  the  time  of  the  acoustic  perturbation  passage  through 
the  layer  with  dimension  of  the  order  of  nk  (the  latter  time 
is  t~nk/vs~  10s).  Thus  the  amplitude  of  the  ion-acoustic 
wave  excited  by  the  acoustic  perturbation  increases  with 
time. 

Above  we  have  determined  some  characteristic  values 
of  the  ion-acoustic  wave.  Generally  speaking,  the  spec¬ 
trum  of  the  ion-acoustic  perturbation  in  k-space  is 
sufficiently  wide.  Indeed,  nk  is  the  length  of  the  acoustic 
pulse.  Therefore  the  acoustic  perturbation  contains 
different  Fourier  components  (not  only  with  k~2nlk). 
For  the  flat  profile  of  the  acoustic  perturbation  in  the 
region  rl<r<r2  (where  r2—ri  =  nk)  and  zero  perturbation 
outside  this  region,  the  Fourier  component  corresponding 
to  the  wave  vector  k  varies  as  (exp(  —  inkk)—  l)/k.  The 
main  contribution  to  the  Fourier  component  (determining 
the  tendency  of  its  decrease  for  large  k)  is  made  by  the 
factor  1  Ik.  This  means  that  not  only  the  waves  with 
to~l  s_l  but  also  with  larger  a>~ Art»s  can  appear  in  the 
spectrum  of  the  ion-acoustic  oscillations. 

We  now  discuss  the  possibility  of  an  appearance  of  ion- 
acoustic  solitons  in  the  ionospheric  plasma.  In  accordance 
with  the  results  of  Alebastrov  et  al.  (1985)  at  the  altitudes 
of  the  order  of  200  km  the  formation  of  ion  acoustic 
solitons  is  possible  as  a  result  of  explosions  with  the  energy 
release  of  about  several  kt  TNT.  In  our  case,  the  electron 
density  perturbation  caused  by  the  acoustic  perturbation 
is  significant:  n'~0.1.  This  implies  the  possibility  of  the 
appearance  of  strongly  non-linear  ion-acoustic  waves. 
Using  the  formula  describing  a  (small-amplitude)  ion- 
acoustic  soliton  propagating  in  the  collisionless  plasma 
(see,  e.g.  Popel  et  al.,  1995)  we  find  the  soliton  width 
Ax  ~  10- 100  cm  for  the  Debye  length  in  the  ionospheric 
plasma  rDe~  l-10cm. 

Let  us  estimate  the  electron  density  variation  caused  by 
an  ion-acoustic  soliton  for  an  arbitrary  energy  release  due 
to  meteoroid  event.  We  use  the  Boltzmann  distribution 
for  the  electron  density  and  assume  that  the  a0-part  of  the 
energy  of  the  acoustic  perturbation  is  transformed  to  the 
ion-acoustic  wave.  For  the  ionospheric  parameters  at  the 
altitude  of  about  200km:  ne~3  x  105cm~3, 
F„~1.5x  103K;  the  energy  of  the  cosmic  body  Ek»  1  kt 
TNT,  oco«3  x  10-4,  we  obtain  1  (the  magnitude  of  a0 
is  chosen  so  that  one  can  obtain  the  magnitude  n'~l 
observed  as  a  result  of  an  explosion  described  by  Ale- 
bastrov  et  al.  (1985)).  Thus  for  the  kinetic  energies  of  the 
Cosmic  body  exceeding  1  kt  TNT  the  variations  of  the 
Electron  density  perturbations  can  be  significant  and  reach 


their  maximum  values  L  Such  strong  perturbations 
can,  in  particular,  influence  radio  wave  propagation. 

Above  we  have  considered  excitation  of  the  electrostatic 
ion-acoustic  oscillations  in  the  ionospheric  plasma 
accompanied  by  a  variation  of  electron  density.  The 
second  important  effect  is  related  to  variations  of  geo¬ 
magnetic  field  strength.  These  variations  can,  in  particu¬ 
lar,  be  created  by  Alfven  waves  excited  as  a  result  of  strong 
actions  on  the  Earth’s  surface  or  on  the  atmosphere  (see, 
e.g.  Galperin  et  al.,  1985).  The  significant  part  of  the 
energy  of  the  Alfven  wave  is  contained  in  its  magnetic 
field.  Significant  magnetic  field  variations  SB~B0  due  to 
the  Alfven  wave  excitation  are  possible  only  for  cosmic 
bodies  with  very  large  kinetic  energy  (of  the  order  of 
100  Mt  TNT). 

Excitation  in  the  ionospheric  plasma  of  sufficiently 
intense  oscillations  caused  by  meteoroid  events  results  in 
the  possibility  of  resonant  acceleration  of  plasma 
particles.  To  describe  processes  similar  to  those  leading  to 
auroral  events  one  has  to  explain  the  generation  of  fast 
electrons  with  energies  from  hundreds  of  eV  to  tens  of 
keV.  Consequently,  the  problem  of  interest  is  to  determine 
the  rate  of  generation  of  electrons  with  energies  of  the 
order  of  1  keV. 

The  process  of  resonant  acceleration  is  described  by  the 
quasilinear  equation  (see,  e.g.  Tsytovich,  1989).  For  the 
case  of  the  ion-acoustic  waves  the  number  of  resonant 
(with  these  waves)  fast  electrons  (v »  vTe,  where  v  is  the 
electron  speed)  is  negligibly  small.  Therefore,  the  process 
of  the  radiative-resonant  acceleration  (Tsytovich,  1989) 
can  be  important. 

Using  the  quasilinear  equation  we  find  that  the  rate  ens 
of  energy  increase  of  electrons  resonant  with  the  waves  is 


8eks  vTe  W 

—  cccoTc - — 

dt  v  n0Tc 


(ID 


where  co^  is  the  electron  plasma  frequency,  W  is  the  energy 
density  of  the  oscillations  in  the  plasma.  The  rate  of  energy 
increase  of  fast  ( v  » i?Te)  electrons  due  to  the  process  of  the 
radiative-resonant  interactions  takes  the  form  (Tsytovich, 
1989) 


deKS 

dt 


(12) 


The  last  estimate  allows  us  to  determine  the  rate  of 
increase  of  the  number  of  fast  electrons  (let  us  remind  that 
we  are  interested  in  consideration  of  particles  with  the 
energy  ~  1  keV) : 


dn?st  1  defast 
dt  ”°efast  dt 


(13) 


The  estimate  for  the  data  R ~ 200  km,  «o~3  x  105cm-3, 
7;~1.5xl03K,  2Tk  ~  1  kt  TNT,  Oo~0.3  x  10“4),  and 
£fest~  1  keV  gives  dnc*st/dt~  10_2n0s_1.  This  magnitude  is 
sufficiently  large  (we  recall  that  the  time  of  existence  of 
the  ion-acoustic  oscillations  is  several  minutes,  i.e.  most 
(or  at  least  a  significant  part)  of  the  electrons  interacting 
with  the  ion-acoustic  oscillations  can  be  accelerated  up 
to  energies  of  the  order  of  1  keV).  Thus  the  process  of 
electron  acceleration  as  a  result  of  the  interaction  of  a 
cosmic  body  with  the  atmosphere  can  be  important  (it 
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can  result  in  excitation  of  atoms  and  molecules  of  the 
ionosphere)  and  should  be  taken  into  account. 


Effects  caused  by  the  presence  of  dust  particles 

If  a  cosmic  body  is  embedded  in  a  dust  cloud  (as  with 
the  comet  Shoemaker-Levy  9)  and  moves  through  the 
magnetosphere  of  the  Earth  then  this  can  lead  to 
additional  electromagnetic  effects  which  are  caused  by  the 
presence  of  dust  particles.  Interaction  of  the  mag- 
netospheric  plasma  with  particles  of  the  dust  cloud  results 
in  charging  of  the  latter.  The  equilibrium  charge  of  a  dust 
particle  is  established  in  the  process  of  equalization  of 
electron  and  ion  currents  on  this  particle,  and  (for  thermal 
equilibrium  distributions)  it  can  be  determined  from  the 
equation  (Tsytovich  and  Havnes,  1 993) 

~exp(— z)  =  — (r+z)  (14) 

^Ti 

where 


a  is  the  characteristic  size  of  the  dust  particle,  —  e  is  the 
electron  charge,  —  Zde  is  the  dust  particle  charge.  In 
hydrogen  plasma,  for  t  =  1  we  have  z  =  2.5.  If  a~  1  fi m 
and  re~  1  keV  (that  is  inherent  in  the  plasma  layer  and 
the  ring  current  of  the  magnetosphere)  then  Zd  is  of  the 
order  of  106. 

The  presence  of  a  dust  cloud  in  the  Earth’s  mag¬ 
netosphere  can  result  in  the  possibility  of  formation  of  a 
global  current  system.  The  effect  of  its  formation  is  due 
to  the  differential  motion  between  the  negatively  charged 
dust  particles  and  the  corotating  (with  the  Earth)  plasma. 
The  dust  grain  moves  across  the  magnetic  field  whereas 
the  magnetospheric  thermal  ions  must  gyrate  along  the 
magnetic  field  lines  and  hence  are  trapped  in  corotation 
with  the  Earth.  The  relative  motion  of  these  two  charge 
carriers  in  the  azimuthal  direction  will  lead  to  a  current. 
The  dust  current  flowing  across  the  magnetic  field  lines 
will  be  mapped  from  the  dust  cloud  to  the  Earth’s  iono¬ 
sphere  via  a  system  of  field-aligned  currents.  The  situation 
is  analogous  to  that  predicted  by  Ip  and  Prange  (1994)  for 
the  case  of  interaction  of  the  comet  Shoemaker-Levy  9 
with  the  Jovian  magnetosphere  and  ionosphere  and  simi¬ 
lar  to  the  planetary  dust  ring/cloud  current  systems  con¬ 
sidered  previously  by  several  authors  (see,  e.g.  Ip  and 
Mendis,  1983;  Goertz  and  Morfill,  1983;  Mendis  et  al., 
1984).  The  evaluation  of  such  a  current  for  the  total  num¬ 
ber  of  dust  grains  in  a  dust  cloud  (Nx  3  x  1020)  and  the 
radius  of  the  dust  cloud  (109cm),  which  are  close  to  those 
related  to  the  comet  Shoemaker-Levy  9  (Ip  and  Prange, 
1994),  and  for  a  space  body  velocity  of  2  x  106cms~'  at  a 
distance  of  about  10  Earth’s  radii  from  the  Earth’s  surface 
gives  us  the  value  of  the  current  of  the  order  of  10s  A. 
Thus  a  sufficiently  intense  current  appears.  The  magnetic 
field  created  by  this  current  at  the  boundary  of  the  dust 
cloud  is  if  ~4  x  10~4G.  This  magnitude  exceeds  the  mag¬ 
netic  field  in  the  plasma  layer  of  the  magnetosphere 
(B~(  1-2)  x  10_4G).  Thus  in  the  vicinity  of  the  dust  cloud 
the  magnetic  structure  of  the  magnetosphere  can  be  viol¬ 


ated.  Furthermore,  the  width  of  the  current  flow  can  be 
comparable  with  the  size  of  the  dust  cloud  (Ip  and  Prange, 
1994).  In  this  case  the  electrodynamical  coupling  of  cur¬ 
rents  flowing  in  different  directions  is  possible.  Elec¬ 
trodynamical  coupling  of  intense  currents  normally  trig¬ 
gers  characteristic  auroral  processes,  including  parallel 
electric  fields,  energetic  particle  acceleration  along  field 
lines  and  precipitation,  together  with  the  generation  of 
radio  waves.  These  processes  would  lead  to  the  detection 
of  an  enhanced  level  of  auroral  UV  and  radio  emission. 

If  the  dust  particles  have  magnetic  moment  (e.g.  if  they 
are  fragments  of  broken  siderolite)  and  one  can  suppose 
that  they  are  distributed  in  space  randomly  then  particle 
acceleration  based  on  the  mechanism  Fermi  is  possible 
(see,  e.g.  Tsytovich,  1989).  The  rate  of  the  energy  increase 
of  particle  is 

<de/d/>  =  4  evic/Lcfi  (16) 

where  e  is  the  energy  of  the  particle,  v  is  its  velocity,  u  is 
the  velocity  of  dust  grains  (i.e.  the  velocity  of  the  cosmic 
body  embedded  in  the  dust  cloud),  L  is  the  mean  distance 
between  the  dust  grains,  c„  is  the  velocity  of  light.  Using 
the  data  of  Ip  and  Prange  (1994)  related  to  the  comet 
Shoemaker-Levy  9 :  the  number  of  dust  grains 
N  ~  3  x  1020  and  the  size  of  the  dust  cloud  r~  109cm,  we 
find  that  the  mean  distance  between  the  dust  particles 
isL~(3A/47rr3)_1/3~2x  102cm.  For«~2x  106cms~'  we 
find  from  equation  (16)  that  electron  can  acquire  a  velocity 
of  the  order  of  the  velocity  of  light  c0  in  a  time  of  the 
order  of  1  s.  This  means  the  possibility  of  creation  of  fast 
particles  which  can  trigger  auroral  processes  before  the 
effective  interaction  of  the  cosmic  body  with  the  atmo¬ 
sphere  (approximately  1  h  before). 


Summary 

In  summary,  the  analysis  of  the  evolution  of  an  acoustic 
perturbation  formed  in  the  atmosphere  due  to  an  impact 
of  a  meteoroid  with  a  kinetic  energy  of  the  order  of  1  kt 
TNT  is  presented.  If  in  the  ionosphere  one  can  neglect  the 
cylindrical  features  in  the  form  of  the  acoustic  per¬ 
turbation  then  the  characteristics  of  the  acoustic  per¬ 
turbation  in  the  ionospheric  plasma  are  close  to  those  of 
the  perturbation  formed  due  to  a  usual  explosion  with  the 
energy  release  coincident  with  that  in  the  meteoroid  event. 
For  the  estimates  we  have  used  experimental  data 
obtained  for  TNT  explosions  on  the  Earth  (experiments 
MASSA).  In  these  cases  an  increase  of  geomagnetic 
activity  was  detected  by  satellites.  Conservation  of  the 
cylindrical  features  in  the  geometry  of  the  acoustic  per¬ 
turbation  in  the  ionosphere  results  in  a  different  action  on 
the  ionosphere  of  a  meteoroid  impact  and  of  the  usual 
explosion  of  the  same  energy  release. 

We  have  discussed  the  possibility  of  excitation  of  ion- 
acoustic  oscillations  in  the  ionospheric  plasma  by  the 
acoustic  perturbation  and  the  possibility  of  an  appearance 
of  ion-acoustic  solitons  in  the  ionospheric  plasma.  The 
process  of  conversion  of  the  acoustic  perturbation  to  the 
ion-acoustic  mode  occurs  at  the  altitudes  of  150-200  km. 
The  time  at  which  the  acoustic  perturbation  reaches  these 
altitudes  is  of  the  order  of  5xl02s.  The  characteristic 
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frequency  of  the  ion-acoustic  oscillations  is  of  the  order 
of  1  s-1.  However  oscillations  with  higher  frequencies  are 
also  present  in  the  plasma.  We  have  presented  the  dis¬ 
persion  law  of  the  ion-acoustic  oscillations  and  their 
damping  rate.  The  time  of  existence  of  these  oscillations 
in  the  ionospheric  plasma  is  several  hundreds  of  seconds. 
High  values  of  the  electron  density  variation  indicate  the 
possibility  of  an  appearance  of  ion-acoustic  solitons  as  a 
result  of  the  meteoroid  events.  For  the  relative  electron 
density  variation  in  the  soliton  of  about  0.1  its  width  can 
reach  values  of  about  100  cm.  We  have  estimated  the 
relative  electron  density  variation  for  arbitrary  energy 
release  due  to  meteoroid  event.  For  kinetic  energies  of 
cosmic  bodies  exceeding  1  kt  TNT  the  relative  electron 
density  variation  can  reach  its  maximum  values  (of  the 
order  of  unity).  Such  strong  perturbations  can,  in  particu¬ 
lar,  influence  radio  wave  propagation.  We  have  estimated 
variations  of  geomagnetic  field  strength  caused  by  the 
meteoroid  events  and  have  given  some  estimates  of  the 
rates  of  electron  acceleration  as  a  result  of  the  interaction 
of  the  cosmic  body  with  the  Earth’s  atmosphere.  Most  (or 
at  least  a  significant  part)  of  electrons  interacting  with 
the  ion-acoustic  oscillations  can  be  accelerated  up  to  the 
energies  of  the  order  of  1  keV. 

Electromagnetic  effects  in  the  Earth’s  magnetosphere 
caused  by  a  cosmic  body  embedded  in  a  dust  cloud  have 
been  also  considered.  The  process  of  dust  particle  charging 
during  the  meteoroid  passage  has  been  taken  into  account. 
The  charge  gained  by  dust  particles  has  been  found.  The 
presence  of  charged  dust  particles  is  shown  to  result  in  the 
formation  of  a  current  system  and  (for  a  sufficiently  large 
concentration  of  dust  particles)  in  violation  of  the  mag¬ 
netic  structure  of  the  magnetosphere  in  the  vicinity  of 
the  dust  cloud.  The  case  of  dust  particles  with  non-zero 
magnetic  moment  has  also  been  considered.  In  this  case 
the  particle  acceleration  based  on  the  Fermi  mechanism 
can  be  important.  The  possibility  of  auroral  processes 
triggered  in  this  case  before  the  effective  interaction  of  the 
cosmic  body  with  the  atmosphere  has  been  shown. 
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